Extracting joint weak values with local, single-particle measurements 
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Abstract 

Weak measurement is a new technique which allows one to describe the evolution of postselected quantum systems. It 
appears to be useful for resolving a variety of thorny quantum paradoxes, particularly when used to study properties of pairs 
of particles. Unfortunately, such nonlocal or joint observables often prove difficult to measure weakly in practice (for instance, 
in optics - a common testing ground for this technique - strong photon- photon interactions would be needed). Here we derive 
a general, experimentally feasible, method for extracting these values from correlations between single-particle observables. 
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In a ground-breaking work, Aharonov, Albert, and Vaidman (AAV) [1] introduced the concept of a "weak mea- 
surement." In contrast to the usual approach to quantum measurment, weak measurements can provide information 
without greatly disturbing the evolution of a system. They are therefore of special interest in the context of post- 
selected subensembles. In particular, they have been used to "resolve" [2] Hardy's Paradox [3]; to interpret and 
extend novel cavity-QED experiments [4, 5] ; to reconcile diverging views on which-path experiments and uncertainty 
[6, 7, 8]; and to address the long-standing controversy over tunneling times [9]. Often, it is of interest to discuss the 
weak value of a nonlocal observable (especially, of course, in proposals to study locality [2, 3, 10]); unfortunately, it is 
usually difficult or impossible to carry out such measurements [2, 11]. In this Letter, we demonstrate that the weak 
values of joint observables may be readily extracted from experimentally realistic measurements. 

In the usual quantum mechanical measurement process, correlations are generated between the shift of the mea- 
surement pointer and the quantum system through strong coupling interactions; this entanglement leads to effective 
collapse. In AAV's weak measurement scheme, the coupling between the measurement device and the system is 
extremely small and their entanglement, and hence collapse, is avoided. The drawback is that on any given event the 
resolution is insufficient to distinguish the different eigenvalues of the observable. Nevertheless, by performing the 
experiment many times on identically prepared systems one can reduce the measurement uncertainty to any desired 
precision. 

AAV considered the case where the weak measurement takes place between two strong measurements called pre- and 
post-selection. Pre-selection of a quantum state is commonplace in quantum mechanics experiments, and amounts to 
the preparation of an initial quantum state, \i) . In a similar manner, post-selection projects out some final quantum 
state, |/). Only upon sucessful pre- and post-selection is the state of the measurement apparatus considered. The 
weak measurement of an obervable A can take place via the same interaction Hamiltonian as a strong (von Neumann) 
measurement, 

H!=gAP x , (1) 

where A is the hermitian operator corresponding to the observable A of the quantum system, g is a (real) coupling 

constant, and P x is the measurement pointer momentum operator. This interaction shifts the mean pointer position 

X^j by an amount Ax = K i^A^ = gT {A^j, where T is the length of time over which the interaction is "on" . The 

main result of AAV's seminal work is that for sufficiently weak coupling strength K, the mean shift in the pointer 

position for a pre- and post-selected subcnscmble is given by Ax = K Re ( A > , where (A) is the expression 

\ / w \ iw 

derived by AAV and termed by them the "weak value" [1]: 




The concept of weak measurement has been controversial since its inception [12]. In contrast to normal quantum 
mechanical expectation values, weak values are not constrained to lie within the extremes of the eigenvalue spectrum 
or to be real numbers. The controversy is only in the interpretation of these strange results, since the predictions 
of the formalism follow directly from quantum mechanics. To date, the predictions of weak measurements have been 
tested only in optical experiments. Ritchie et al. performed the first such test in an optical Stern-Gcrlach experiment 
[13], based on some theoretical proposals [1, 14]. More recently, we have performed an experimental realization of 
the Quantum Box Problem [15, 16]. In these experiments, the single-particle observables of interest were coupled to 
the transverse displacement of a light beam; this transverse position acted as the measurement pointer. 

Optical experiments are very limited in the kinds of interactions that are available to experimentalists at the 
quantum level. One can realize an effective interaction of the form of Eq. 1, for single-photon polarization [13] or 
optical path measurements [16]. However, there are more complicated types of observables for which it is difficult 
to find a strong interaction Hamiltonian of the desired form. An example of such an observable would be the joint 
probability of one particle being found in region A and a second particle being found in region B. Using AAV's 
approach, one would require a Kerr-like nonlinearity at the single-photon level to couple the value of this joint 
property into a shift of the measurement pointer. Such nonlinear effects are very difficult to come by in optics 
[17, 18, 19, 20], and thus optical experiments might seem suited to only the simplest weak measurement problems. 
More complicated theoretical proposals, such as the ones involving Hardy's Paradox [2, 3], nonlocality of single particle 
[10], and extensions of the Quantum Box Problem [21] are extremely difficult to implement if they require optical 
nonlinearities at the quantum level, or other multi-particle or nonlocal interactions. One theoretical proposal has been 
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made for carrying out such joint weak measurements with trapped ions, relying on the natural, strong, interparticle 
interaction [11]; to date, no such quantitites have been measured in the laboratory. 

In this work, we show that it is experimentally possible to measure these complicated joint weak values using 
local, single-particle interaction Hamiltonians by studying the correlations between two different pointer variables 
(either two different one-dimensional pointers or a single two-dimensional pointer). We derive an explicit formula 
for finding the weak value of an observable which is the product of two operators; this is sufficient to test all of the 
aforementioned theoretical proposals in optical systems. (The theoretical approach is general enough that it can be 
extended to products of N operators by using N pointer dimensions.) 

We begin by rederiving AAV's result for the deflection of a measurement pointer through the weak measurement 
procedure using the interaction Hamiltonian from Eq. 1 . Throughout our calculations wc assume that the operators 
A and P x commute, which is a very natural assumption since one operator corresponds to a property of the quantum 
system and the other to the measurement apparatus. 

The quantum system is initially prepared in the state We consider a measurement pointer initially prepared 
in a gaussian state centered at zero: 



1 



27TCT 



exp 



x 
"4^2 



(3) 



where a is the rms width of the probability distribution. We now calculate the expected displacement of the 
measurement pointer when the quantum state is post-selected in the final state |/), using the rule 



X) = 

fi 



P(f)X 

P(f) 



(4) 



where \XJ is the expectation value of X given a successful post-selection of the final state |/) , and P(f) is the 

projector |/) (/| onto the final state [9]. This amounts to finding the expectation value of the operator O x = \ f) (f \ X 
and renormalizing. We employ the Heisenberg equation of motion to propagate the expectation value of O x forward 
in time to first order: 



O x {t)) ~ (0.(0)) + 



H,O x 



(5) 



This first-order approximation is only valid for weak values of the coupling K. Our measurement pointer was originally 
centred at zero, and therefore the first term on the right-hand side vanishes. The expectation value of O x to first-order 
is 



O x (t) 
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The values of the pointer integrals are ( (j> 
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d in truncating 



= -ih/2 and 

the time-dependence of the system at first order only if the higher-order terms are small. Second-order terms will 
contain integrals with quadratic terms in the momentum and linear in position and are therefore zero by symmetry; 
however third-order terms will be proportional to l/er 2 . These terms will contribute negligibly to the expectation 
value of the weak value provided Ka -C a, where a is the characteristic scale of A; that is to say that the uncertainty 
in the pointer position is much larger than the strength of the interaction - this was AAV's condition for weakness. 
We now renormalize to obtain the conditional expectation value of X given pre-selection and postselection, 
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Thus we see that the shift of the pointer, conditioned on appropriate postselection is indeed proportional to the weak 

value (A) from Eq. 2: 
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It is well known that weak values are not constrained to be real numbers. The imaginary part of the weak value 
yields a shift in the variable conjugate to the pointer position (i.e. P x ), related to the back-action of the measurement 
[9]. Using the same initial state of the measurement apparatus, one can go through a similar calculation for the 
expectation value of the operator O p = \f) (f\ P x . The result of this calculation for the imaginary part of the weak 
value is: 



Im 
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r- ■ (10) 

w n K 

Note that, unlike the real part, the imaginary part of the weak value depends quadratically on the width of the pointer 
distribution, a. This means that for more and more ideal weak measurements (i.e. as a — > oo) the actual physical 

shift in the pointer momentum, \P X ) , will tend to zero. 



We now proceed to calculate an expression for joint weak values. Instead of employing a nonlinear optical Hamil- 
tonian to measure a joint weak value (^AB^ of two observables A and B (with corresponding hermitian operators 

A and B), we use the local single-particle interaction Hamiltonian 

H = gAAP x +g B BP V7 (11) 

where g A and gs are coupling constants and P x and P y are the momentum operators conjugate to two commuting 
pointer operators X and Y. Our system is prepared in the quantum state \i) and the measurement apparatus in a 
state |$) = \<j) x ) \(j) y ) described by a two-dimensional gaussian wavepacket, 
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where a x y are the rms widths of the probability distribution in x and y, respectively. 

In this case, we are interested in the joint weak value i^-Bj , and therefore need an expression for a pointer 

deflection that contains contributions from deflections in the x-direction and in the y-direction after the successful 
post-selection of the quantum system into the state |/). By analogy with the previous calculation, we find the 
expectation value of the operator, O xy = \f) (f\XY. We use the Heisenberg equation of motion recursively to 
expand the right-hand side to second order: 



O xy (t)) ~ (O xy (0)) + 
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The zeroth-order term, (O xy (0)\ , in the expectation value is zero since |$) is even with respect to both x and y, 



(14) 



while O is odd in both. The first-order commutator is given by: 



H, xy ] = (g A AP x + geBPy) \f) (f\ XY - \f) (f\ XY (g A AP x + g B BP y ) . 



Each term in this expression is even with respect to either x or y but not both, and therefore has a vanishing 
expectation value for the pointer state |5>). Therefore, to first order, since the measurement apparatus can receive a 
kick in either the x-direction or the y-direction but never both, there is no shift in the expectation value of O. We 
require the second-order term, where kicks in both x and y can appear, to find a nonzero value. The second-order 
commutator contains 16 terms of which only 8 contain the product of P x and P y and can therefore have nonzero 
expectation values. These remaining 8 terms can be simplified by collecting 2 identical pairs of terms to 
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BAP x P y O xy + ABP x P y O xy 1AP x O xy BP y 
c 2BP y O xy AP x + O xy ABP x P y + O xy B AP x P y 



We further simplify this result by evaluating the integrals for the measurement apparatus, 
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where K x , y = gA,Bt in analogy with our one-dimensional definition. We again employ the notation (xY^ to mean 

the expectation value for XY given a successful post-selection. This conditional expectation value requires division 
by the success probability for normalization such that 
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where and B\y are the weak values of the individual operators as defined by AAV (Eq. 2), and can be extracted 
from experimental measurement as shown in the previous theory section. We again are making an approximation 
when truncating our evolution at second order. Third order terms are identically zero due to symmetry; however 
fourth-order terms are not. Similar to the single weak values, these higher order terms contribute negligibly provided 
that K XyV are sufficiently small with respect to the widths <J x , y , in direct analog with the condition we imposed in the 
one-dimensional case. We rearrange Eq. 19 to obtain an expression for the weak value of (AB + BA)/2, 
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where the subscript W refers to the weak value. If A and B are commuting observables (as must be the case, for 
instance, if they operate on disjoint subsystems), then AB is Hermitian and the left-hand side is equal to Be(AB)w, 
the desired result: 



Re (AB" 



' K X Ky 



Re (vC(*)J- 



(21) 



For AB non-Hcrmitian, the AAV coupling Hamiltonian 7i = gABP would be unphysical, and only the symmetric 
form of Eq. 20 could correspond to an observable (assuming A and B Hermitian and g real). Note that other operator 
orderings have also been studied in certain contexts [22]. 

Through an analogous calculation we can obtain an expression for the imaginary part of the weak value of the 
operator AB by calculating the first nonzero term in the expectation value of the operator O xpy = \f) (f\ XP y under 
the same interaction Hamiltonian from Eq. 11. The imaginary part of the weak value is: 



Ab + bA\ 
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(22) 



which can again be further simplified on the left-hand side if A and B commute. As we noted in the single weak value 



case, the physical implication (i.e. the shift in (^XP y ) ) of the imaginary part of the weak value becomes smaller and 



21 and 22 are experimentally measurable in the same experiment. The 
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w 
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smaller in the limit as a — > oo. 

Note that all of the quantities in Eqs 

quantities (^XY^ and (^XP y ^ arc extracted from the two-dimensional pointer deflections, and 

can be measured using the method from Eq. 9 and Eq. 10. 

Since AAV's original work on weak measurement, their predictions have been demonstrated experimentally in 
several optical experiments. In a number of circumstances, it has become clear that weak measurements of nonlocal 
and/or multi-particle observables should be of great interest. The interaction Hamiltonians required to accomplish 
such measurements using the original AAV approach are often experimentally inaccessible, e.g., demanding strong 
interactions between individual photons. We have developed an alternative method for extracting the desired joint 
weak values that can be implemented using only local, single-particle Hamiltonians. In addition to being necessary for 
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the resolution of certain quantum paradoxes [2, 10, 16], the ability to discuss the correlations present in postselected 
systems should prove indispensible for investigating the role of postselection in quantum information [23]. The new 
approach described here will be essential for the broad experimental application of joint weak measurements. 

We would like to thank Jeff Lundeen and Morgan Mitchell for helpful discussions. This work was funded by the 
U.S. Air Force Office of Scientific Research (F49620-01-1-0468), Photonics Research Ontario, and NSERC. 
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